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Abstract 

The Lagrange mesh method is a very accurate and simple procedure to compute eigenvalues and 
eigenfunctions of nonrelativistic and semirelativistic Hamiltonians. We show here that it can be 
used successfully to solve the equations of both the relativistic flux tube model and the rotating 



Oh| string model, in the symmetric case. Verifications of the convergence of the method are given. 
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I. INTRODUCTION 



The Lagrange mesh method is a very accurate and simpleprocedure to compute eigenval- 
ues and eigenf unctions of a two-body Schrodinger equation [l|, |2j, y| . The trial eigenstates are 
developed in a basis of well-chosen functions, the Lagrange functions, and the Hamiltonian 
matrix elements are obtained with a Gauss quadrature. This method can be extended to 
treat very accurately three-body problems, in nuclear or atomic physics jm . Recently, it has 
also been successfully applied to a two-body spinless Salpeter equation p . The idea of this 
work is to adapt the Lagrange mesh method to solve the complicated equations of both the 
relativistic flux tube and the rotating string models. 

The relativistic flux tube (RFT) is a phenomenological model describing the mesons. 
It relies on the assumption that the quark and the antiquark are connected by a straight 
color flux tube carrying both energy and momentum. The quarks are considered as spinless 
particles in the original version of the model 0, Q, Isj]. The RFT reproduces the linear 
Regge trajectories, and reduces to the usual Schrodinger equation with a linear confinement 
potential in the nonrelativistic limit. We will consider here the particular case of mesons 
composed of two equal quark masses. The equations of motion of the symmetric RFT model 
are given by two coupled nonlinear equations: one defining the Hamiltonian and the other 
defining the orbital angular momentum. These equations depend on a quark transverse 
velocity operator and their solutions will be obtained by the use of an iterative procedure 
similar to the one proposed in Ref. 

The rotating string model (RS) also describes the mesons. It is derived from the QCD 
Lagrangian and is characterized by the fact that it contains auxiliary fields. The equations 
of motion for this model are similar to the equations of motion of the RFT model jjj . 
In the symmetric case, it has been showed that the RS is classically equivalent to the RFT 



if the auxiliary fields are correctly eliminated 



asymmetric case , provides a clear physical interpretation for the characteristic variables 
of the RS model. 

The Lagrange mesh method is explained in Sec. |TIJ In Sec. the relativistic flux tube 
and the rotating string models are described. Then, it is shown, in Sec. II VI how the Lagrange 
mesh method can be applied to solve the equations of motion of these models. After some 
remarks, given in Sec. about the numerical and physical parameters, the results are 



ll|. This result, extended recently to the 
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presented in Sec. IVII and the reliability of our numerical method is checked. Finally, some 
concluding remarks are given in Sec. IVHI 



II. LAGRANGE MESH METHOD 

A Lagrange mesh is formed on iV mesh points Xi associated with an orthonormal set of 
indefinitely derivable functions fj(x) ODD A Lagrange function fj(x) vanishes at all 
mesh points but one; it satisfies the Lagrange conditions 

.!)(■'■;) = A, ; %, (1) 

The mesh points x^, the zeros of a particular polynomial, and the Aj are connected with a 
gauss quadrature formula 

/ g(x) dx « ^2 X k g(x k ) , (2) 
Ja fc=i 

used to compute all the integrals over the interval [a, b}. 

As we consider only radial equations, this interval is [0, oof, leading to a Gauss-Laguerre 
quadrature. The Gauss formula (j2J) is exact when g(x) is a polynomial of degree 2N — 1 
at most, multiplied by exp(— x). The Lagrange-Laguerre mesh is then based on the zeros 
of the Laguerre polynomial L^{x) of degree N [J. An explicit form can be derived for the 
corresponding regularized Lagrange functions 

fi(x) = (-iyx; 1/2 x{x - x^L^x) e~ x ' 2 . (3) 

To show how these elements can be applied to a physical problem, let us consider a 
Hamiltonian H = T(p 2 ) +V(r), where T(p 2 ) is the kinetic term and V(r) a radial potential 
(h = c = 1). The calculations are performed with trial states \ip) given by 

N 

fc=l 

where 

(r \h) = fj ^ Y Ur). (5) 

I is the orbital angular momentum quantum number and the coefficients Ck are linear 
variational parameters, h is the scale parameter chosen to adjust the mesh to the domain 
of physical interest. We define r = hx, with x a dimensionless variable. 



We have now to compute the Hamiltonian matrix elements. Using the properties of the 
Lagrange functions and the Gauss quadrature (J2J), the potential matrix is diagonal. Its 
elements are 

(MVirm^Vihx^, (6) 

and only involve the value of the potential at the mesh points. As the matrix elements are 
computed only approximately, the variational character of the method cannot be guaranteed. 
But the accuracy of the method is preserved 

The kinetic energy operator is only a function of p 2 . Let us define the corresponding 
matrix, 

P* = (W 2 \fj)- (7) 

It is shown in Ref. ^] that, using the Gauss quadrature and the properties of the Lagrange 
functions, one obtains 



where 



(-1) 1 'fax,) V^Xi + x^Xi-Xj) (i^j), 

Prij = \ ( 9 ) 

(12x?)- 1 [4+(4JV + 2)z i -a;?] (i = j). 

Now, the kinetic energy matrix T(P 2 ) can be computed with the following method jlj: 

1. Diagonalization of the matrix P 2 . If D 2 is the corresponding diagonal matrix, we have 

P 2 = SD 2 S~\ (10) 
where S is the transformation matrix. 

2. Computation of T(D 2 ) by taking the function T of all diagonal elements of D 2 . 

3. Determination of the matrix elements = (fi\T(P 2 )\fj) in the Lagrange basis by 
using the transformation matrix S 

T(P 2 ) = ST(D 2 )S-\ (11) 

This procedure can easily be generalized to the case of an arbitrary function F of any given 
matrix M, in order to compute F(M) (provided the calculation is relevant). Note that such 
a calculation is not exact because the number of Lagrange functions is finite. However, it 
has already given good results in the semirelativistic case, where T(p 2 ) = \/p 2 + m 2 



The eigenvalue equation reduces to a system of N mesh equations 



A? 



[Ty + V(hxi)5ij — E5ij] Cj = with Cj = ^h\ju{hxj) 



(12) 



where u(r) is the regularized radial wave function. The coefficients Cj provide the values 
of the radial wave function at mesh points. But contrary to some other mesh methods, the 
wave function is also known everywhere thanks to Eq. (J1J. 

III. THE MODELS 

A. The relativistic flux tube 

In the original RFT model jfj , the meson is composed by two spinless particles - a quark 
and an antiquark - which move being attached with a flux tube. This tube is assumed to 
be linear with a uniform constant energy density a and carries angular momentum. A tube 
element has only a transverse velocity. The system rotates in a plane with a constant angular 
velocity u around the center of mass, assumed to be stationary. If is the distance between 
the ith quark and the center of mass, and if we define fj = dri/dt the radial velocity oh the 
ith quark, then the quark speed is given by vf = r\ + vf ± , where Vi± = tori. We also assume 
that the energy density of the extremities of the flux tube is modified of a negative constant 
C/2, in order to take into account possible boundary effects due to the contact between the 
tube and the quark. Further, we consider that the quarks can interact via V(r) taking into 
account a short-range potential (a one-gluon-exchange process, for instance). These two 
extra terms are discussed in Ref. j8|. The Lagrangian C of the meson is given by 



d = -mai 1 -a dr'^ 1 - -7^, (14) 

where mj is the constituent mass of the zth quark, ji — (1 — vf )~ l l 2 and ji± = (1 — vf^)^ 1 ^ 2 . 

In the following, we will only consider the symmetric case, mi = 1712 = m. Then, n = r2, 
r = 2ri, t>ij_ = t> 2 j_ = v±. The corresponding quantized equations of the system are 0,0] 




(13) 





(15) 



r 




(16) 
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where I is the orbital angular momentum, {A, B} = AB + BA, 4x 2 f(x) = arcsinx — 
x\J\ — x 2 , W r = \/p 2 + m2 ; an d p 2 . = —\j^zr. The operator v± commutes neither with 
r nor with p r operators. These equations reduce to a spinless Salpeter equation with the 
potential ar + V(r) + C when I = 0, and to a Schrodinger equation with the same potential 
in the nonrelativistic limit. The general case (mi ^ 7712) is detailed in Ref. Q. 



B. The rotating string 

Starting from the QCD Lagrangian and writing the gauge invariant qq Green function 
for confined spinless quarks in the Feynman- Schwinger representation, one can arrive at the 
Nambu-Goto Lagrangian, which describes two quarks with masses mi and m 2 , attached by 
a string of energy density a. With the straight line ansatz and the introduction of auxiliary 
fields (einbein fields) to get rid of the square roots appearing in this Lagrangian, one can 
obtain the Hamiltonian 



* = 1 

2 



Pr +rn\ , Pr + m 2 2 2 2 f 1 dp ''' 



+ — + Hx + /x 2 + a V / — + / + r 

Ml ^2 Jo V JO a 3T 2 



+ V(r), (17) 



where 



a 3 = /ii(l - C) 2 + M2C 2 + /' d/3 (/3 - C) 2 ^ (18) 

The potential V(r) takes into account interactions not simulated by the rotating string. We 
do not consider here a contribution coming from a constant potential C, as in the RFT model. 



L = y/ t{t + 1) and £ defines the position i? M of the center of mass: i? M = (£1^ + (1 — Qx2ni 
where is the coordinate of the ith quark. The auxiliary fields \L\ and \Li can be seen as 
effective masses of the quarks, while the auxiliary field v can be interpreted as an effective 
energy density for the string. 

We are interested here in the resolution of the symmetrical case. When m\ = m 2 = m, 
then C = 1/2 and ^1 — ^2 — A*- Defining 

y = T L -> (is) 

2a 3 r 

one can eliminate v by a variation of the Hamiltonian. This extremal field u Q reads 

z/ = ar =. (20) 

Jl -4y 2 ((3- 1/2)2 
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By replacing v by v Q in the Hamiltonian (JT7j) and the relation (|19p. we obtain the following 
equations for the symmetrical rotating string [lO 1 ] 

+ ^arcsiny - yV 1 - V J , (21) 



ar 2 ar ' Ay 2 

nrp 1 _|_ jyi^ Qnp 

H = — h/iH arcsinw + ay 2 + V(r). (22) 

» y 



It has been shown in Ref. pjj that the extremal value of p giving 8H/8p = is 



-=J%$- (23) 



Moreover, the replacement of p by po in Eqs. f)21jl and ({22} gives exactly the symmetrical 
RFT equations (|T3j) and (|T5|l. with y equal to The RS model with all its auxiliary fields 
eliminated is thus equivalent to the RFT model in the classical symmetrical case. This is 
also true when (mi ^ mo), as shown in Ref. Il2j. 

nn 

Here, we use the RS model with the auxiliary field p not eliminated, as in Refs. !<). 10. In 
these papers, the parameter p is considered as a real parameter and not as an operator. But, 
to avoid eventual singularities in the value of this auxiliary field when y is classically close 
to 1, we introduce explicitly the dependance of p in y, through the following substitution 

" "> -7=> ( 24 ) 

v 1 -v 

where p is a real number. Such an expression is inspired by the result (j2HJ). The quantized 
equations of the symmetrical rotating string are thus 

= P^T=^ + 2 {rJ{y)h (25) 

1 f o 9 / ^1 1 + y 2 a f arcsinyl Tr , , , , 

# = - (p? + m 2 , Vl - y 2 } + P^f= 2 + 2 \ r ' "V/ + (r) ' ( } 

where 4x 2 /(x) = arcsinx — x\f\ — x 2 like for the RFT model. 

A particular solution depends on the value of this parameter p. Following Refs. SQ, the 
physical value of p minimizes the mass of the state. The mean value (p) — (p/ a/1 — y 2 ) can 
be considered as a constituent mass for the quark, depending on the state. These equations 



reduce to a Schrodinger-like equation with the potential ar + V(r) when I = JjJ, and to 
a true Schrodinger equation with the same potential in the nonrelativistic limit. 
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IV. RESOLUTION 



A. The relativistic flux tube 



The main purpose of our work is the resolution of the symmetrical flux tube equations 
(fT3Jl and (JTfij) using the Lagrange mesh method. To do this, we have to compute the matrix 
elements of the different operators in the Lagrange basis. As we consider a radial problem, 
we will use a Gauss-Laguerre quadrature. So, the corresponding Lagrange functions will be 
given by Eq. Q. Let us define the different matrix elements we need to know 



As 



arcsm v 



B, 



(fMfj) 



Av 



fj 



(Mi±\fi) 



Gij = (fi\v±7±\fj) , Sij 



V* = (W{r)\fj) • 



With these notations, Eqs. (fTHjl and (fTHjl read 

A = {G, D} + a{B, F} + CG, 



fi 



31 i 



arcsm v± 

V I 



(27) 



(28) 



H = {T,D} + -{B,S} + CT 
where we have used the approximative closure relation, 

N 



V, 



i. 



k=l 



(29) 



(30) 



to compute a product of two matrices. 

The matrix elements Aij, Bij, and Vij are easy to compute, thanks to Eq. (jHJ). Moreover, 
Eq. © gives us an analytical expression for p^ij, from which we can deduce the matrix 
elements by using the procedure described in Sec. [TTJ The same procedure will allow us to 
compute Fij, Gij, Sij and once the matrix elements of v± are known. The determination 
of these matrix elements can be achieved by an iterative process, described here: 



1. Equation ()28j) can be rewritten as 
G 

where 



-{P^- 1 } - -{CD' 1 } - -DGD- 1 - -D~ l GD, 
2 l J 2 1 J 2 2 



P = A-a{B,F}. 



(31) 



(32) 



This equation is symmetrized to ensure that G is hermitian. It is worth noting that 
P = P(G) since F = F(G). Starting from an known matrix G k at the kih step, P k 
can be computed and we obtain a new matrix G k ' with Eq. (|31j). 

2. This iterative process would diverge if we choose G k+1 = G k> . So, we introduce a new 
parameter e < 1 and define G k+1 = eG k ' + (1 — e)G k . 

3. At each step k, the N eigenvalues {^if of the operator v± are computed. The 
iteration procedure ends when 

1 N 

i=l 

where rj is a fixed tolerance. 

Once we have reached the convergence for G, we are able to compute S and T, which are 
now seen as functions of the matrix G rather than the matrix elements of the operator v±. 
The Hamiltonian can then be computed and diagonalized. 

Actually, the final matrix G is practically independent of the initial one G°. However, 
the faster way to reach the convergence is to develop Eq. (fT3j) at the first order in v± and 
to choose the matrix G given by this development. At the first order, v±^± ~ v±, and 

G° « VWTT) Q{fi, D} + ^ + 1 1) 1 . (34) 
Let us note that a relevant starting matrix is obtained even if m = 0. 



(fe+l) (k) 



(k+1) 



< V, 



(33) 



B. Rotating string 



1. Lagrange mesh method 



The resolution of the RS with the Lagrange mesh method is similar to that of the RFT. 
Indeed, using the previous definitions (|27)) with y instead of v±, and defining 

' ■-■<- 



Qv = (fi | ^A 3 F| ./• ) , Yii = (ft ^ /,) 



/,), E ij = {f i \p l r +m 2 \f j ) 



Eqs. fl23|) and (J2Sj) are given by 



G = —(A — a{B,F}) , 
2p 

H = ^{E,Q} + pY + ^{B,S} + V. 



(35) 

(36) 
(37) 



Like for the RFT, we need to compute the matrix of the operator y to completely know the 
Hamiltonian. We will do this by an iterative process on G given directly by Eq. ()36|) . with 
an initial value, obtained after a first order development, given by 

G° = vW+T) (pB + ^B 2 ) _1 . (38) 

The last step in the resolution of the RS equations is always to find the value of the real 
number p realizing the minimum mass of a particular state. This extremal value is different 
for each state. 



2. WKB method 

Contrary to the case of the RFT, the operator p 2 appears only in the equation defining 
the Hamiltonian for the RS. This makes possible a solution of Eqs. f!25|) and (ffilj) by a WKB 
method. 

First, let us examine the case I = 0. The RS equations reduce then to a spinless Salpeter 
equation of the form (p = p since y = 0) 

H=t—- + p + ar + V(r), (39) 

where 

P 2 = Pr + ^2 " ( 40 ) 

In the WKB method, L = £ + 1/2. Consequently L 2 = 1/4 here, and we obtain 

p 2 . = pM — p 2 — m 2 — par — — pV(r), (41) 

where M is the meson mass. We have then to compute r + and r_ the two physical zeros of 
the classical quantity p 2 . Finally, the resolution of the Bohr-Sommerfeld condition 

Prdr = tt > ( 42 ) 

followed by a minimization of M with respect to the parameter p gives the mass of the state 
whose quantum numbers are I and n. 
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When C ^ 0, the WKB formulation of the classical RS equations (J2Tj) and with the 
substitution pijl. reads 



+ TT arcsin 2/ -yv 1 -^ 2 )> ( 43 ) 



ar 2 ar a/1 - y 2 % 2 

L / 1 - \~ xji^ civ 

M = ~(pt + m 2 ) J I - y 2 + p , H arcsinw + V(r). (44) 

p v ; y 

The first one implicitly defines a function y = y(r, £, p), which can be numerically computed. 
We can then formally write 

p clv 1 -\- y 2 
pz = — (M — V(r) — arcsiny) — p 2 — m 2 . (45) 

^i^W y i-y 2 

The rest of the resolution is now identical to the previous case i = 0. 



V. SET OF PARAMETERS 
A. Physical parameters 

In this paper, we are mainly interested in the capacity of our method to give accurate 
solutions of the coupled equations for both RFT and RS models. But, in order to com- 
pare our results with previous studies and to use our method with physical parameters in 
interesting ranges, we will use the values of physical quantities from the models la and Ic 
developed in Ref. (see Table HJ). Both models possess a coulomb term with three values 
of the strength, depending on the quark content of the meson: k^i for heavy-light system, 
Khh for heavy-heavy system, and Ku for light-light system (light quark: u, d, s; heavy quark: 
c, b). 



B. The scale parameter 

The Lagrange mesh method provides us a direct picture of the wave function at the mesh 
points. The best results are thus obtained when the mesh covers the main part of the wave 
function and the last mesh point is located in the asymptotic tail. That is why we are 
interested in an adequate determination of the scale parameter h. Since the method is not 
variational, no extremum of the mass can be expected for a defined value of h. A good value 
for this quantity is given by h = r a /xN, where is iVth zero of the Laguerre polynomial 

11 



(the last point of the mesh), and r a represents a distance for which the asymptotic tail of 
the wave function is well defined. If a; at is well known, r a is not. We show here how such a 
quantity can be estimated. 

A typical evolution of the computed masses for different values of h is presented in Fig. 
The existence of plateaus shows that the method does not require the knowledge of precise 
values of the scale parameter. A simple estimation will be sufficient, even to obtain accurate 
results. 

For given quantum numbers, a system of two massless quarks is expected to have the 
maximal spatial extension, and so it could give an upper bound of the parameter h. First, 
we analyze the problem for the RS equations when I = 0. These equations reduce then to 
a spinless Salpeter Hamiltonian, which reads 

H A = — + p + ar. (46) 
P 

We fix V(r) = 0, since the asymptotic behavior is controlled by the confinement. The 
solutions have the following analytical forms (n = 0, 1, . . . ) ^Jj 

P 



Eno(p) = (-) (s n )+p, (47) 



Uno{r) = [pa) /0 S (48) 



where Ai(s) is the Airy function and s n its nth zero, given by the approximate formula 



|U1 



37T . 3 X 



-| 2/3 

(49) 



or 



Replacing p by its extremal value p„o, 

Pno = Va [^^J ' ( 50 ) 
we have 

u n0 (r) Ai ^/a (^~Y^j r + s^j . (51) 

When s ~ 5, Ai(s) is about 0.02% of its maximal value. Consequently, a good estimation 
of r a is given by 

_ \ 1/4 

S " ^ r a + Sn = 5. (52) 
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At this point, we are able to compute a "physical" value for h when £ = 0. The extension 
of the wave function increases with the angular momentum. The simplest way to simulate 
such an increase is to compute r a with the relation 



This crude estimation of h is satisfactory because it is always located in the plateau. More- 
over, we will use it in both RFT and RS methods, because of the classical equivalence 
between these two theories. 

C. Numerical parameters 

The accuracy of the solutions depends mainly on two parameters: the number N of mesh 
points (basis states) and the value of the tolerance r\ on the eigenvalues of the operator v±. 
For instance, a relative error on meson masses around 10~ 5 can be reached with N > 30 and 
7] < 10~ 6 . The accuracy can be increased by using greater values of N and smaller values 

of Tj. 

If the value of the mixing parameters e is too high, the iterative process diverges. The 
best value of e is chosen as the largest value for which the process converges. It depends 
on the case considered, as shown in Table |Hj It clearly appears that the iterative process 
does not converge easily with the RFT equations, especially when the quarks are massless. 
About 700 iterations are needed in this case, and 400 when mj \fa > 1. However, the RS 
solutions converge faster, and one can reach the convergence after about only 40 iterations. 

VI. RESULTS 

A. Relativistic flux tube 

We have computed with the Lagrange mesh method the solutions of the RFT equations 
for models la and Ic from Ref. js| (see Table HJ). All the masses are computed with N = 30, 
1] = 10 -6 , the scale parameter h is estimated thanks to Eq. (|53jl. and the parameter e 
is taken from Table |H] Meson masses are presented in Table IIHI with the corresponding 
ones computed with the method developed in Ref. [8|, relying on a harmonic oscillator 




(53) 
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basis. Experimental data are given in order to show that the parameters used are physically 
relevant. 

The results of both methods are compatible. Nevertheless, the masses computed with the 
Lagrange mesh method are always smaller than the masses computed with the harmonic 
oscillator method, although the Lagrange mesh method is not variational. Our method 
provides thus a better convergence of the results. The improvement is especially important 
for light quark masses. Differences between the two methods vanish when the quark mass 
increases. 

It is worth noting that the masses computed with method of Ref. j^J are strongly de- 
pendent of the values chosen for the oscillator length. So a supplementary minimization on 
this parameter, for each state, is necessary to obtain the optimal value of a mass. This is 
not necessary with the Lagrange mesh method since it is nearly independent of the scale 
parameter (see Fig. QJ). 

The small differences between the masses obtained with the Lagrange mesh method and 
the harmonic oscillator method are a strong indication that our method works well. But we 
want another test. It will be given by the study of the Rotating string model. 

B. Rotating string 

Solutions of the RS equations computed with the Lagrange mesh method (numerical 
parameters as in Sec. IVI A|) and the WKB approximation are presented in Table IIVI The 
masses are obtained using the set la of parameters (see Table |TJ), for a pure string without 
coulomb-like potential. 

The two methods to solve the RS equations lead to very close results. This shows that 
the semiclassical approximation is efficient in this case, but also that the Lagrange mesh 
method works correctly. Fig. El shows the existence of a minimal mass for a particular value 
Pq of the parameter p in the RS equations. In our calculations, p has been determined to 
the nearest 10 MeV, and is the same in the two methods with that precision. An accuracy 
below 1 MeV is then reached for the masses. 
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C. Comparison between RFT and RS 



If the RS model is classically equivalent to the RFT model once the auxiliary fields are 
correctly eliminated, the two models should not give the same results when a real parameter 
p is kept in the RS equations. In a previous study 11], some results have been obtained 
about the equivalence between a spinless Salpeter Hamiltonian Hss and a corresponding 
Hamiltonian with auxiliary field Ha'- the eigenvalues of Ha are upper bound of the eigen- 
values of Hss 13 with relative differences around 7% for the lowest states. We know that 
the RFT and RS equations reduce respectively to eigenvalues equations for Hamiltonians 
Hss and Ha for a vanishing angular momentum. It should be interesting to see if there is 
the same kind of relation between the masses for the RFT and RS models when i ^ 0. 

Another result can be expected: once we know an eigenfunction \ipRs) of the RS Hamilto- 



nian for the extremal value p , we are able to compute the effective mass p = p (1/ a/1 — y 2 ) 
for this state. This quantity should be approximately equal to the mean value Prft — 
(V(.Pr + m2 ) 7(1 ~~ v ±) / f° r t ne corresponding state \iPrft), due to the equivalence between 
the two models via Eq. 

Our results are given in Table IVl The masses for both RFT and RS models are computed 
with the Lagrange mesh method for the same parameters as in Sec. IVI Bl The RS masses 
are always upper bound of the RFT masses with relative differences around by 7%, as in 
the limiting case of vanishing angular momentum. We also see that Prft ~ Ho as expected. 
We can finally notice that the results of the two models are closer and closer when the mass 
of the constituent quark increases, because the RFT and the RS model posses and common 
nonrelativistic limit: the Schrodinger equation with a linear potential. 



VII. CONCLUSIONS 

We have shown in this paper that the Lagrange mesh method solves successfully the 
equations of the relativistic flux tube model in the symmetrical case. The masses obtained 
are in good agreement with a previous resolution in a harmonic oscillator basis [8] . But the 
Lagrange mesh method is more efficient, due to its independence of the scale parameter used 
to fit the size of the trial states. Moreover, a better convergence is reached. This proves the 
validity of our method. 
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We have also solved the equations of the symmetrical rotating string model with the 
Lagrange mesh method and with the WKB approximation. The masses computed with 
these two procedures are very close, showing that the Lagrange mesh method correctly 
works, and that the WKB approximation is efficient here. If we compare the masses given 
by the relativistic flux tube and the rotating string models, we find relative differences around 
7% for the lowest states, as expected because the two models are classically equivalent. This 
point is a last confirmation of the relevance of the Lagrange mesh method to solve the 
relativistic flux tube equations. 
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TABLE I: Two sets of physical parameters for the RFT and the RS models, from Ref. [8J (n = u 
or d). 





la 


Ic 


m n (GeV) 





0.233 


m s (GeV) 


0.317 


0.416 


m c (GeV) 


1.456 


1.658 


a (GeV 2 ) 


0.151 


0.169 


C (GeV) 


o 




/>■// 


1.016 


0.539 


«w 


0.698 


0.467 




0.544 


0.500 



TABLE II: Approximate optimal values for the parameter e in different cases. 







e 






RFT 




RS 


ps 


0.005 




0.1 


>1 


0.01 




0.1 
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TABLE III: Meson masses for the RFT model, with two sets la and Ic of parameters from Ref. [8j, 
computed using the Lagrange mesh method (Lag.) and a previous technique relying on an harmonic 
oscillator basis (HO) 8]. The experimental masses (Exp.) are given, without error, for information. 



Mass (GeV) 



(n + l) 25+1 Lj 


Exp. 


HO (la) 


Lag. (Ia) 


HO (Ic) 


Lag. (Ic) 


nn 1 3 S\ 


0.771 


0.781 


0.762 


0.774 


0.773 


1 3 P 2 


1.318 


1.310 


1.300 


1.320 


1.319 


1 3 D 3 


1.691 


1.654 


1.643 


1.689 


1.676 


2 3 S 1 


1.465 


1.450 


1.415 


1.427 


1.424 


2 3 Po 


1.810 


1.841 


1.832 


1.797 


1.794 


ss l 3 Si 


1.019 


0.988 


0.968 


1.010 


1.010 


1 3 P 2 


1.525 


1.540 


1.534 


1.517 


1.515 


l 3 D x 


1.854 


1.881 


1.877 


1.867 


1.865 


2 3 Si 


1.680 


1.671 


1.641 


1.644 


1.641 


2*P 2 


2.011 


2.053 


2.047 


1.994 


1.991 


cc l 3 Si 


3.097 


3.131 


3.130 


3.116 


3.115 


1 3 P 2 


3.556 


3.528 


3.527 


3.542 


3.542 


1 3 D 3 


3.770 


3.788 


3.788 


3.820 


3.820 


2 3 Si 


3.686 


3.666 


3.663 


3.664 


3.661 


2 3 D X 


4.159 


4.128 


4.128 


4.165 


4.164 
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TABLE IV: Meson masses for the RS model computed with the Lagrange mesh method and the 
WKB approximation. The interaction la from Ref. is used, but without the Coulomb potential. 
The extremal value po, to the nearest 10 MeV, of the parameter p is given in both cases. 



Lagrange mesh WKB 







Mass (GeV) 


Po (GeV) 


Mass (GeV) 


Po (GeV) 


nn 




1.289 


0.32 


1.294 


0.32 




1 3 P 2 


1.581 


0.16 


1.589 


0.16 






1.960 


0.49 


1.963 


0.49 


cc 


1 3 S 1 


3.492 


1.55 


3.493 


1.55 






3.731 


1.52 


3.730 


1.51 






3.916 


1.65 


3.917 


1.65 
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TABLE V: Meson masses computed with the Lagrange mesh method within the RFT and the RS 

n . 

models. The interaction la from Ref. |8( is used, but without the Coulomb potential. The values 



of the effective masses 


Ho and ijlrft are , 


given. 










Relativistic flux tube 


Rotating 


' string 






Mass (GeV) 


Hrft (GeV) 


Mass (GeV) 


u n (GeV) 


nn 


l 3 Sl 


1.228 


0.308 


1.289 


0.32 




1 3 P 2 


1.543 


0.323 


1.581 


0.29 




1 3 D 3 


1.825 


0.342 


1.860 


0.32 




2 3 5i 


1.832 


0.460 


1.960 


0.49 




2 3 P 2 


2.071 


0.498 


2.155 


0.49 


ss 


l 3 5x 


1.507 


0.486 


1.536 


0.49 




1 3 P 2 


1.809 


0.523 


1.838 


0.52 




1 3 ^3 


2.078 


0.593 


2.103 


0.56 




2 3 5i 


2.070 


0.612 


2.142 


0.62 




2 3 P 2 


2.294 


0.647 


2.343 


0.64 


cc 


l 8 Sl 


3.486 


1.555 


3.492 


1.55 




1 3 P 2 


3.723 


1.594 


3.731 


1.58 




1 3 D 3 


3.931 


1.625 


3.937 


1.62 




2 3 5i 


3.902 


1.631 


3.916 


1.65 




2 3 P 2 


4.081 


1.661 


4.094 


1.66 
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FIG. 1: Typical evolution of meson masses for the RFT model with the scale parameter h: IS, IP 
and 2S states for the isospin 1 mesons, computed with the parameters la from Ref. ^. Formula (|5M|) 
gives h = 0.21 GeV -1 for the IS state; this value is correctly located in the plateau. 



FIG. 2: Meson masses for the RS model, computed with the Lagrange mesh method, versus the 
parameter p: IS, IP and ID states for the ss system, computed with the interaction la from 
Ref. 8], but without the Coulomb potential. 
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h(GeV ') 




(GeV) 



